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ABSTRACT 

Tarnowsky, Terence J. Ph.D., Purdue University, May 2008. Long- Range Multiplic- 
ity Correlations in Relativistic Heavy Ion Collisions as a Signal for Dense Partonic 
Matter. Major Professor: Rolf P. Scharenberg. 

A dense form of matter is formed in relativistic heavy ion collisions. The con- 
stituent degrees of freedom in this dense matter are currently unknown. Long-range, 
forward-backward multiplicity correlations (LRC) are expected to arise due to mul- 
tiple partonic interactions. Model independent and dependent arguments suggest 
that such correlations are due to multiple partonic interactions. These correlations 
are predicted in the context of the Dual Parton Model (DPM). The DPM describes 
soft partonic processes and hadronization. This model indicates that the underlying 
mechanism creating these long-range multiplicity correlations in the bulk matter is 
due to multiple partonic interactions. 

In this thesis, long-range multiplicity correlations have been studied in heavy ion 
(Au+Au and Cu-|-Cu) and hadron-hadron (pp) collisions. The behavior has been 
studied as a function of pseudorapidity gap {Arj) about rj = 0, the centrality, atomic 
number, and incident energy dependence of the colliding particles. Strong, long-range 
correlations {Arj > 1.0) as a function of Arj are found for central collisions of heavy 
ions at an energy of ^/sNN = 200 GeV. This indicates substantial amounts of dense 
partonic matter are formed in central heavy ion collisions at an energy of ^/SNN = 
200 GeV. 



1. Introduction 

1.0.1 Quark Searches 

Searches for quarks in high-energy accelerator and cosmic ray experiments, and 
even experiments similar to the Millikan oil-drop experiment, have not observed free 
quarks. However, in deep inelastic scattering events, where a high-energy electron 
scatters off a proton or neutron, the cross section of wide-angle scattering was incon- 
sistent with an electron interacting with a point-like proton. In the 1960s, high-energy 
electron experiments became available to probe the fundamental structure of nucleons 
via deep inelastic scattering. At larger momentum transfers (Q^), it is possible to 
probe short-range features. The physical existence of quarks and gluons inside pro- 
tons and neutrons was discovered in these high-energy inelastic electron scattering 
experiments. 

1.0.2 Experimental Basis for Quantum Chromodynamics 

The basic tenets of the strong interaction and its characterization by quantum 
chromodynamics (QCD) have been based on experimental results. The existence of 
hadrons consisting of three quarks in the same spin state suggested that quarks must 
possess an additional distinguishing characteristic. One example, the A"''^ (or A^) 
resonance consists of three u (d) quarks in the same spin state (|). Since u and 
d quarks are fermions, they will obey the Pauli Exclusion principle. As such, the 
ground state wave function must be antisymmetric with respect to the interchange of 
any two quarks in coordinate space. However, based on the known properties of the 
A++ (or A^) it appeared that the wavefunction was symmetric. The problem was 
solved by the addition of a new quantum number, which forms a fully antisymmetric 



Name 


Symbol 


Mass {GeV/c'^) 


Quark Content 


Pion 


7r+ 


0.140 


ud 


Kaon 


K+ 


0.494 


us 


Phi 


$ 


1.019 


ss 


J/^ 


J/^ 


3.097 


cc 


Proton 


P 


0.938 


uud 


Neutron 


n 


0.940 


udd 


Delta 


A++ 


1.232 


uuu 


Lambda 


A" 


1.116 


uds 


Xi 


"ep 


1.315 


uss 


Omega 


Q- 


1.672 


sss 



Table 1.1 
A listing of some hadrons, their masses, and quark content [?] . 



wave function. The concept was proposed by Oscar Greenberg, who discussed quarks 
in terms of a new quantum number. It was Murray Gell-Mann who named the 
new quantity "color". In the case of A++ (or A^), the three quarks have different 
colors and the overall wave function of the resonance is antisymmetric with respect to 
interchange of any quark. The quarks combine in such a fashion to make all baryons 
and mesons color neutral objects (Table fTTTI) . 



In the highest energy scattering events quarks were never observed to exist out- 
side of their parent hadrons. This indicated that the force between two quarks must 
increase as the distance between them increases. The presence of a massless par- 
ticle to mediate the force between two quarks was introduced. This particle was 
referred to as a gluon. The gluon is the analogue to the photon, a massless boson 
that mediates the electromagnetic force between charged objects in quantum elec- 



Flavor 


Symbol 


Current Mass {GeV/c'^) 


Electric Charge 


up 
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~ 0.003 


2 
3 


down 
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~ 0.006 


1 
3 


charm 


c 


1.25 


2 
3 


strange 


s 


~ 0.10 


1 
3 


top 


t 


~ 174 


2 
3 


bottom 


b 


4.20 


1 
3 



Table 1.2 
A listing of quarks, their current masses, and electric charge [?]. 



trodynamics (QED). Unlike the photon, which carries no electric charge, the gluon 
cannot be neutral with respect to the color charge. If this were the case, the potential 
energy for a quark-quark or quark- ant iquark state would yield an average potential 
that for mesons would be negative (attractive), while that for baryons would be pos- 
itive (repulsive). If gluons are colored, the potential for both mesons and baryons is 
attractive. Since they are colored, gluons not only interact with quarks, but other 
gluons as well. Their colored nature leads to gluon self-interaction. 



Since color charge is a conserved quantity, and the emission or absorption of a 
gluon leads to a color change, gluons can be considered as consisting of a color — 



color (anti-color) charge. There are nine possible combinations of color — color, 
which leads to the prediction that there should exist nine gluons. However, the 



ninth combination is actually a linear combination of all color — color, and therefore 
colorless. If this gluon existed, colorless baryons would be able to emit these particular 
gluons and interact with each other via a short-range gluon force, which is not the 
case experimentally. 



The potential between two quarks as a function of distance can be parameterized 

by, 

Ca 
V(r) = -^^^ + kr (1.1) 

r 

where as is the strong couphng constant, r the separation between the quarks, and 
k a constant of about 1 GeV/fm. The form of this potential is such that at small 
separation distances (r ^ 0), the linear term vanishes and one is left with the form 
of an attractive Coulomb potential. The strong coupling constant a^ has been ob- 
served to decrease with increasing momentum transfer. This property is referred to 
as "asymptotic freedom". For large values of r, on the order of 10~^^ m (0.1 fm), the 
linear term dominates. As separation distance increases, an infinite amount of energy 
becomes required to separate the two quarks. However, as the energy in the color 
field between the two quarks increases with the separation distance, it will eventually 
exceed some threshold where it becomes more energetically favorable to create a, q — q 
pair from the vacuum than to continue separating the quarks. This process forms 
two mesons, i.e., colorless objects. This property of the strong interaction makes the 
observation of isolated colored objects impossible. A schematic drawing of the chro- 
moelectric flux between a quark-antiquark pair is shown in Figure 11.11 The six known 
quark flavors are shown in Table 11.21 along with their electric charge and estimates 
of their current mass (quark mass when deconfined) [?] . The constituent mass is the 
effective mass of a quark in a binding potential and is larger than the current mass. 

1.1 Phase Diagram of Strongly Interacting Bulk Matter 

The quantum chromodynamic theory of strongly interacting matter predicts a 
phase transition between hadronic matter and quark-gluon degrees of freedom as the 
temperature or density of hadronic matter is increased. A schematic representation 
of the proposed QCD phase diagram is shown in Figure [L2l [?]. 
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Figure 1.1. Flux lines between a quark- ant iquark pair, which collapse into 
a tube between the two particles as their separation distance increases [?]. 
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Figure 1.2. A proposed phase diagram for QCD showing temperature as 
a function of baryon chemical potential (/is) [?]. High-energy heavy ion 
collisions probe the regime of low fiB and high temperature. The small 
curve denotes the nuclear liquid-gas phase transition [?,?,?]. 
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Normal nuclear matter exists at low temperature and a baryon chemical potential 
of approximately 922 MeV, corresponding to an energy density of ~ 0.15 GeV/fm^. 
At non-zero temperatures in this region the dominant contribution to the vacuum 
pressure is from pions, the lightest mesons. In a pion gas, there are only three 
hadronic degrees of freedom. Neglecting the pion mass, the pressure of an ideal pion 
gas can be expressed as [?], 

P ■ = —T^ (1 2) 

J- pion gas qn \^-'^J 

which is the pressure of a Stefan-Boltzmann blackbody with three degrees of freedom. 
Here it is assumed that the relation between pressure and energy density is P = |e. 
The number of degrees of freedom in a state of matter consisting of deconfined quarks 
and gluons is much greater than that of a pion gas. In a system of deconfined light 
quarks (u, d, and s) and gluons there are a total of N ?« 48 partonic degrees of freedom, 
with 16 bosonic (Nfeoson) degrees of freedom and 36 fermionic (Nfermion) degrees of 
freedom. The bosonic degrees of freedom arise from the gluons that have 2 spin states 
and 8 color combinations (2*8). The fermionic degrees of freedom are contributed 
by both the quarks and antiquarks (2) with 2 spin states, 3 color states, and n^ flavor 
states each (2 * 2 * 3 * Uf). The factor of | in Equation ll.3al arises from differences 
in Bose-Einstein (gluons of spin 1) and Fermi-Dirac (quarks of spin ±|) statistics in 
the Boltzmann factor, nj = 3 if only u, d, and s quarks are considered. In an ideal 
quark-gluon plasma (QGP) with massless quarks, the system can be modeled as a 
massless, relativistic gas with energy density e = 3P. 

7 

-'' -'■^boson ' '^-'■^ fermion yi.Oa) 

O 

21 
N = lQ + —nf (1.3b) 

Pqgp = ^T' (1.3c) 

When the pressures of the pion gas and the QGP are equal, both phases should 
exist in equilibrium. Calculations have been carried out using lattice QCD methods 
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Figure 1.3. Lattice QCD calculations for pressure as a function of tem- 
perature. At ^ = 1 (~ 170 MeV, the critical temperature) there is an 
increase in pressure, potentially indicating a phase transition [?]. 



to determine the deconfinement temperature, Tc. Lattice QCD models quarks and 
gluons on a discrete lattice. From these results, Tc is predicted to lie at ~ 170 ± 10 
MeV, corresponding to temperatures on the order of 10^^ Kelvin. 

Figures 11.31 and 11.41 show lattice QCD calculations for the pressure and energy 
density normalized by T^ as a function of temperature [?]. The calculations are for 
three light quark flavors, two light and one heavy quark flavor, 2 light quark flavors, 

T_ 

and the pure gauge case, i.e., with infinite quark masses. The factor e '"s , taking 
into account the mass of the strange quark irts., accounts for the difference between 
the two light flavor and the "2+1" case. 

At approximately T = T^ (170 MeV) there is a rapid change in the energy density 
of the system. For T > 1.5Tc, the energy density (related to the number of degrees of 
freedom) appears to be nearly constant, but below the Stefan-Boltzmann limit for an 
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Figure 1.4. Lattice QCD calculations for energy density as a function of 
temperature. At ^ = 1 (~ 170 MeV) there is a rapid change in energy 
density, indicating a rapid increase in the number of degrees of freedom 
and a potential phase transition [?]. The arrows indicate the estimated 
temperatures in collisions at SPS, RHIC, and LHC 
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ideal gas of quarks and gluons. Since the pressure rises more slowly than the energy 
density it has been suggested the QGP that may be produced at RHIC (up to ~ 2 
Tc), may still consist of clusters of strongly interacting quarks and gluons [?]. Thus, 
at temperatures greater than 2Tc the coupling may continue to decrease, eventually 
(T>>Tc) yielding a weakly interacting plasma. 

Lattice QCD calculations with realistic quark masses indicate the presence of a 
critical point (region) [?]. The exact location of the critical point is dependent on 
the quark masses used in the calculation. The critical point is shown in Figure 11.21 as 
the point lying on the phase boundary between hadronic matter and a QGP, at fiB 
(baryon chemical potential) equal to normal nuclear matter density. It is expected 
that for a transition to the left of the critical point, there exists a smooth crossover 
region, whereas to the right of the critical point a first order phase transition is 
expected. In the crossover regime a second order phase transition is expected, where 
there is no discontinuity in the entropy, i.e., no latent heat. For the quark-gluon 
to hadronic phase transition, it is possible that chiral symmetry is restored in the 
partonic phase. If chiral symmetry is exact, the pion would be massless. The mass 
of the pion (140 MeV) is non-zero because chiral symmetry is not exact and leads to 
non-zero quark masses {rriu^d 7^ 0). 

1.2 Signatures of the Quark Gluon Plasma 

There are several suggested signatures of quark gluon plasma formation: J/\E' 
suppression [?], strangeness enhancement [?], direct photons, Hanbury Brown-Twiss 
(HBT) determined source radii, etc. [?,?]. Two other observables that can indicate 
QGP formation are related to the bulk, collective nature of the system and the pres- 
ence of energy densities high enough to support a deconfined quark-gluon system. 
These two measurements involve transverse elliptic fiow of the bulk particles and the 
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suppression of high pt particles in the medium. However, none of these signatures 
are conclusive evidence for the formation of quark-gluon matter. 

1.2.1 Strangeness Enhancement 

In pp collisions the production of strange particles is suppressed compared to 
particles containing light {u and d) quarks. Equal yields of n, d, and s quarks are not 
measured. One explanation for the suppression is the higher production threshold 
due to the mass of the s and s quark pair, if thermal production goes as e ~ , where 
iTLq is the mass of the produced quark. Strangeness can be enhanced due to the higher 
temperature of a quark-gluon medium, which leads to a lower effective quark mass 
(related to chiral symmetry restoration) [?]. A deconfined system of quarks and gluons 
can also lead to the production of s and s quark pairs by gluon fusion. Strangeness 
enhancement has been seen experimentally at AGS, SPS, and RHIC, across a wide 
range of energies in both nucleus-nucleus and proton-nucleus collisions. This makes 
the signature of strangeness enhancement an ambiguous signal of the transition to 
quark-gluon matter. 

1.2.2 High-pT Particle Suppression 

The production of a dense, colored medium in heavy ion collisions was predicted 
to lead to the attenuation of high-p^ partons that traverse the medium. If a high den- 
sity, colored medium is produced, a colored parton that passes through the medium 
will interact with the medium constituents via gluon Bremsstrahlung [?]. The final 
consequence of this is a suppression of high px {> 1 GeV) particles. This effect has 
been seen in Au+Au collisions at an energy of ^/JNN = 200 GeV, but is absent in 
d+Au at the same energy [?]. The energy loss of a parton via gluon radiation through 
a particular medium is dependent on the transport coefficient q [?]. The calculation of 
q yields the same energy loss for an ideal pion gas and an ideal quark-gluon plasma, 
at equal energy density [?, ?]. Independently, jet quenching cannot distinguish a 
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deconfined system of quarks and gluons from a hadronic system with similar gluon 
content [?]. 

In the study of nucleus-nucleus collisions there may be initial state nuclear effects 
that influence hadron production. This effect can be quantified using the nuclear 
modification factor. The nuclear modification factor is generally defined as [?], 

Rab{pt) = 7^^-, / , ,2 (l-^) 

where Tab = {Nbin) l^^neias accounts for the collision centrality. This is the ratio 
of inclusive charged hadron yields from A+B collisions to pp, normalized by the 
number of binary nucleon-nucleon collisions. In the absence of nuclear effects, if 
particle production scales with the number of binary collisions (A^bm) then Rab = 1- 
Figure 11.51 shows Rab for central and minimum bias d+Au and central Au+Au data 
at an energy of ^Jsnn = 200 GeV. However, there is a large suppression in particle 
production at high-p^ in central Au+Au data. 

1.2.3 Open Questions 

While the signatures discussed in Section 11.21 provide evidence that dense matter 
is formed in heavy ion collisions, there remain a number of fundamental questions 
that must be addressed. The signals of deconfinement and the bulk properties of 
the system (thermalization, number of degrees of freedom, etc.) remain important 
open questions. Models including partonic or hadronic systems have been used to 
fit the data. The analysis of forward-backward (FB) multiplicity correlations has the 
potential to distinguish between hadron-hadron and parton-parton interactions. Both 
model independent and model dependent arguments suggest that long-range forward- 
backward correlations are due to multiple partonic interactions. This measurement 
can provide a clear signal that partonic degrees of freedom are involved. 
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d+Au FTPC-Au 0-20% 
d+Au Minimum Bias 




8 10 

Pj (GeV/c) 



Figure 1.5. Nuclear modification factor (Rab) for central and minbias 
d+Au and central Au+Au collisions at an energy of ^/snn = 200 GeV, 
showing the strong suppression of high-p^ particle production in central 
Au+Au collisions [?]. 



13 
1.3 Forward-Backward Multiplicity Correlations 

Correlations that are produced across a wide range in rapidity are thought to 
reflect the earliest stages of a heavy ion collision, free from final state effects 



The study of forward-backward (FB) multiplicity correlations has a long history. 
Studies were carried out in e+ — e^, fi^ — p, p — p, p — p, and other experiments 
[?,?,?,?,?,?,?, [2n]- Correlations that extend over a long range in pseudorapidity {i]) 
have the potential to probe the early stages of heavy ion collisions. Forward-backward 
correlations have been characterized by the forward-backward correlation strength, 
b, the slope extracted from a linear relationship between the average multiplicity 
measured in the backward rapidity hemisphere (< Ni, >) and the multiplicity in the 
forward rapidity hemisphere, Nj. This relationship was predicted theoretically, seen 
in hadron-hadron experiments, and expressed as [?,?], 

< Nb{Nf) >=a + bNf (1.5) 

In this definition, the correlation strength h can be positive or negative with 
a range of |6| < 1. This maximum (minimum) represents total correlation (anti- 
correlation) of the produced particles separated in rapidity. 6 = is the limiting 
case of entirely uncorrelated particle production. Experimentally, the slope of h in 
hadron-hadron experiments is found to be positive. The intercept of Equation 11.51 a 
is related to the number of uncorrelated particles. 

The correlation strength can be expressed as the ratio of the covariance of the 
forward-backward multiplicity and the variance of the forward multiplicity. This is 
done by performing a linear regression of Equation 11.51 and minimizing R^, 

n 

^' = E[^^. -(« + &%)]' (1-6) 
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The minimization condition implies -^ = -|^ = 0. This produces a set of equa- 



tions that can be expressed in matrix form, 



-1 



n 



ENl~ EN 



U 



1=1 \i=i 

i=l 1=1 1=1 

When expanded to solve for b, the matrix in Equation 11.71 yields. 



T.Nb.T.Nl-T.NnT.Nf.N,^ 



i=l i=l 4=1 i=l 

n n n 



;i-7) 



n 



n 



n 



nY.NfM.-llNuY.N,. 

i=l i=l i=l 



n 
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the covariance of A^;, and N ^ normalized by the variance of N f. Equation 11.81 can 

be expressed in terms of the following calculable average values, 
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with [12 



< NfN, >-<Nf><N,>= dy dy' [< N{y)N{y') >-< N{y) >< N{y') >] 

Jy<0 Jy'>0 

(1.10a) 
<Nf>-<Nf >2= f dy f dy' [< N{y)N{y') >-< N{y) >< N{y') >] 

Jy>0 Jy'>0 

(1.10b) 



1.3.1 Short-Range Correlations 



Short-range correlations are correlations that extend over a small range of pseu- 
dorapidity {\ri\ < 1.0). Short-range correlations are due to various short-range order 
effects [12] . These effects can include particles produced from cluster decay, resonance 
decay, or jet correlations. The particles produced in a single inelastic collision are 
known to only exhibit short-range correlations [?]. 
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1.3.2 Long- Range Correlations 

Long-range correlations are correlations that extend over a wide range in pseudo- 
rapidity, beyond |?7| > 1.0. The presence of long-range correlations is a violation of 
short-range order. Short-range order is expected to hold as long as "unitarity con- 
straints are neglected" [12]. In the approximation of short-range order, only single 
scattering is considered. Therefore, quantum mechanical probability is not conserved, 
since it is possible to have multiple scattering terms. 

The presence of short-range order in each inelastic collision produces short-range 
correlations, discussed in Section 11.3.11 The consideration of unitarity leads to the 
existence of multiple inelastic scattering, in addition to the single scattering that 
determines the short-range correlations. Multiple inelastic elementary scatterings 
are the source of the long-range correlation [?]. Due to the possibility of multiparton 
scatterings in hadron-nucleus and nucleus- nucleus collisions compared to pp, if a long- 
range correlation exists at a given energy it should be enhanced in collisions involving 
nuclei [?]. Long-range correlations have been seen in high-energy (a/s/vw > 1 TeV) 
pp collisions where multiparton excitations are seen in high multiplicity events [?,30j. 
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2. Experiment 
2.1 RHIC 

The Relativistic Heavy Ion Collider (RHIC) complex (Figure 12. ip at Brookhaven 
National Laboratory (BNL) consists of a succession of coupled accelerators that pro- 
vide RHIC with fully ionized heavy ions, or polarized protons. RHIC is a supercon- 
ducting hadron collider with two concentric rings (colloquially referred to as "Blue" 
(clockwise ion revolution) and "Yellow" (counter-clockwise ion revolution), 3.8 km in 
circumference. As Figure 12.11 shows, the rings are not circular, but consist of several 
arc sections mated to straight, insertion sections. The center of the insertion sections 
contains the intersection points. There are six intersection points around the RHIC 
ring. Four are (or were) occupied by experiments, BRAHMS at 2 o'clock, STAR 
at 6 o'clock, PHENIX at 8 o'clock, and PHOBOS at 10 o'clock. As of 2007, only 
PHENIX and STAR are actively taking data. The superconducting magnets in each 
ring operate at a magnetic rigidity, Bp = 8400 Tesla ■ m, which yields a kinetic en- 
ergy of 100 GeV/A for Au ions, with an ^ ratio of 2.5. The magnetic field strength 
is approximately 3.5 T. At this rigidity, lighter ion species can reach 125 GeV/A, 
while protons, with -^ = 1, can be accelerated to 250 GeV. Due to the presence of 
two independent rings, RHIC has the capability of colliding particles ranging from 
protons {A = 1) to gold (Au) {A = 197) in both symmetric and asymmetric collisions 
(e.g., deuteron+Au). Collisions involving yet heavier particle species up to uranium 
(U) will be possible with the upgrade of an Electron Beam Ion Source (EBIS). 

The injection of heavy nuclei into RHIC requires several successive steps. For 
Au, negatively charged ions are created from a pulsed sputter ion source, partially 
stripped of electrons with a foil at the high voltage terminal of the Tandem Van de 
Graaff, and accelerated by the second stage of the tandem to an energy of 1 MeV/A. 
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Figure 2.1. The Relativistic Heavy Ion Collider (RHIC) complex [T]. 
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After leaving the Tandem and passing through yet another stripping foil, the ions 
have a charge of +32. They are then transferred to the AGS Booster synchrotron, 
where their energy is increased to 95 MeV/A. Upon exit from the Booster the ions 
are further stripped to a charge of +77, leaving only two electrons of the initial 79 in 
the gold atom. These are injected into the Alternating Gradient Synchrotron (AGS), 
the final boosting stage prior to injection into RHIC and accelerated to an energy of 
8.86 GeV/A. They then leave the AGS and are stripped of their final electrons as 
they travel through the transfer line to RHIC, and are then fully ionized to a charge 
of +79. Once in RHIC, the ions are accelerated to their final colliding energy and 
stored for several hours. 



2.2 STAR 

The Solenoidal Tracker at RHIC (STAR), is one of the two large experiments at 
RHIC (Figure 122] PJ)- STAR is a large acceptance detector, covering a pseudorapidity 
range of \rj\ < 1.8, with additional coverage at forward pseudorapidity |?7| ~ 4.0. The 
main tracking detectors in STAR have full azimuthal symmetry, covering A0 = 27r. 
The large rj and acceptance, symmetry about r^ = 0, and accurate tracking makes 
STAR ideally suited for characterizing the bulk particles {pr < 2.0 GeV) and for 
performing detailed correlation studies. 

2.2.1 Time Projection Chamber 

The main detector at STAR is the Time Projection Chamber (TPC), seen in 
Figure l2T3l The TPC is a large cylinder, 4.2 meters in length by 2 meters radius. The 
inner radius of the TPC is 50 cm from the beam axis. Using a TPC, 3-dimensional 
reconstruction of charged particle tracks from a collision is possible. The TPC allows 
charged particle tracking, momentum determination, and identification from energy 
loss (dE/dx) over a range of ± 1.8 units in pseudorapidity {rj, see Appendix A for 
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Figure 2.2. The STAR detector. 



definitions of kinematic variables). Tfie momentum range of particles identifiable by 
the TPC extends from 0.1-1 GeV/c, while the total range of measured momentum 
extends from 0.1-30 GeV/c. The TPC gas volume is filled with PIO, a mixture of 90% 
argon (Ar) and 10% methane {CH4). A central cathode membrane at a potential of 
28 kV, with the end caps at ground potential, supplies the drift electric field in the 
TPC. Electrons liberated from gas atoms by the passage of charged particles through 
the TPC drift longitudinally toward the end caps, where the tracking information is 
read out. The drift time is approximately 40 /xsec, corresponding to a drift velocity 
of ~ 5.5 cm//isec. The minimum two-track resolution of the TPC is 2 cm. The TPC 
(along with all STAR tracking detectors) lies inside the field generated by the STAR 
magnet. At its maximum, this uniform, longitudinal magnetic field has a magnitude 
of 0.5 T. 
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A schematic diagram of the STAR Time Projection Chamber. [2] 
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Figure 2.4. Cross section of the STAR detector and subsystems. 



2.2.2 Forward Time Projection Chambers 



Two additional TPCs he at forward rapidities covering the phase space 2.5 < 
\ri\ < 4.0 over the fuU azimuth. These Forward Time Projection Chambers (FTPCs) 
are unhke the main TPC in that they utihze a radial drift electric field, compared 
to the longitudinal drift for the main TPC. This was a necessary result of the small 
extent of the FTPCs, which are only ?^ 1 m in length. The radial drift improves 
two-track resolution to ~ 2 mm, an order of magnitude better than the main TPC. 
The drift electric field is created between a cylindrical inner cathode and the outer 
wall at ground potential. The electrons ionized by the passage of charged particles 
through the active gas volume drift to the outer radius of the FTPCs where their 
signals are read out by 9600 pads. The maximum radial drift distance is 22.32 cm. 
The FTPC gas volume consists of a 50-50% mixture of Ar and carbon dioxide (CO2). 
Further information regarding the work done for the FTPCs is explained in Appendix 
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2.2.3 Trigger and Data Acquisition 

The trigger detectors at STAR consist of the central trigger barrel (CTB), the 
zero-degree calorimeters (ZDCs), and the Beam-Beam Counters (BBCs). The CTB 
measures charged particles within 77 = ± 1 with full azimuthal coverage. The ZDCs 
are located up and downstream of the STAR detector, at the position of the RHIC DX 
magnets, the point where the incoming (outgoing) ion beams are brought together 
(split). The ZDCs measure spectator neutrons from the collision that propagate with 
the beam momentum. The location of the primary collision vertex can be determined 
from time of flight differences in a coincidence signal in the two ZDCs. The BBCs 
provide an additional coincidence measure for triggering and are also used to monitor 
beam luminosity. The ZDCs and BBCs provide complementary vertex determination 
and are used in the majority of minimum bias triggers. 

The data acquisition system (DAQ) acquires data at rates up to 100 Hz. Due to 
the large multiplicities in heavy-ion collisions, raw data file sizes can be as large as 
200 MB per event. At typical trigger rates, this is far too large for DAQ to handle. To 
reduce the file size, zero- suppression is applied for each detector. Events are written 
to tape storage at the RHIC Computing Facility (RCF) at a rate of 30-50 MB/s. 
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3. Theoretical String Model 

3.1 Color Strings/Chains 

Multiparticle production at high energies can be characterized by a string model 
(from a chromoelectric flux tube) via Schwinger pair production, which was first de- 
rived for e^e~ collisions [3J. It was later generalized to the case of hadron collisions, 
including heavy ion collisions Pl[5]. Between a produced q — q pair is a constant, 
chromoelectric field that manifests as a tube (Figure 11.11) . In this field additional 
q — q pairs can be generated. The creation of new q — q pairs screens the original 
field. This leads to pairing of new quarks with the original ones. The process can be 
repeated many times, leading to multiple hadron production of mesons and baryons. 

In heavy ion (A + A) collisions, it is expected the number density of these strings is 
much higher than in e~^e~ or pp collisions. In pp or p + A collisions, there may only be 
one string formed. In A + A collisions there can be many strings. Immediately after 
the collision of two heavy ions a chromoelectric field is formed between the receding 
nuclei. These fields are confined in a small transverse area, forming color flux tubes 
(also referred to as "strings" or "chains"). In central (full overlap) collisions, some 
strings may overlap with each other. 

3.2 Dual Parton Model 

The Dual Parton Model (DPM) [6] (and the similar Quark Gluon String Model ^) 
is intended to describe soft hadronic physics in nucleus-nucleus, hadron-nucleus, and 
hadron-hadron collisions. Unlike hard processes, which can be treated perturbatively, 
soft physics occurs in the QCD regime of strong coupling. The DPM has been used to 
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Figure 3.1. Particle emission from string decay. Each q — q pair weakens 
the color field [5]. 
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Figure 3.2. Dominant two string particle production mechanism in the 
DPM for pp interactions at high energy [8]. 



predict several observables of high-energy collisions including: multiplicity, pseudo- 
rapidity (r/), and p^ distributions; particle ratios; charge distributions; KNO scaling; 
heavy particle production; nuclear stopping; J/\l' suppression; strangeness enhance- 
ment; and long and short-range multiplicity correlations in pseudorapidity {rf). 

In pp collisions, the dominant production mechanism is a two string process. 
Figure 13.21 shows this process. The two strings are comprised of the valence quarks 
and therefore create strings with a quark and diquark at the ends. The strings are 
formed by soft gluon exchange between partons. A one string process dominates in 
cases where a q or g in the projectile can annihilate with its antiparticle in the target, 
such as in pp interactions. Higher order contributions also exist. The second order 
contribution is from a four string process, shown in Figure 13.31 The additional strings 
are formed between sea quarks and antiquarks. 
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Figure 3.3. Four string particle production mechanism in the DPM for 
pp interactions [8]. This includes the two string process (Figure [32]) with 
the addition of two strings of the type q — q as shown. 
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3.2.1 Hadronization in the DPM 

Particle production in the DPM is due to the hadronization of chains (or strings) 
stretched between partons. Similar to the Schwinger particle production mechanism, 
(where at some threshold separation distance between a. q — q pair, a new q — q pair 
is created from the vacuum), the chains in the DPM hadronize at some threshold to 
form q ~ q pairs. In high-energy particle collisions, the number of chains produced 
is equal to twice the number of inelastic nucleon-nucleon collisions. The main as- 
sumption in the particle production process put forth by the DPM is that the chains 
hadronize independently. Though the superposition of chains is possible, this inde- 
pendent hadronization scheme is the basis for all DPM predictions. In pp, the single- 
particle inclusive cross section is described by the superposition of two chains [6] . For 
any one cut Pomeron, particles are produced with only short-range rapidity correla- 
tions. It is the superposition of several cut Pomerons (appropriately weighted) that 
contribute to the fluctuation in the number of strings and produces long-range rapid- 
ity correlations. The Pomeron was introduced as the particle exchanged to explain 
the increasing scattering cross-section in high energy experiments at energies greater 
than ^ 20 GeV [SHU]. The exchange of other particles, such as mesons, leads to 
a predicted decrease in the cross-section as a function of energy, whereas Pomeron 
exchange predicts an increasing cross-section with energy. The Pomeron is predicted 
to consist of quarks and gluons combined in such a way that the particle carries no 
color charge. 

3.2.2 Correlations in the DPM 

From the Dual Parton Model, rapidity correlations are described by the average 
multiplicity in a backward rapidity interval, < A^b >, as a function of total multi- 
plicity in a forward rapidity interval, Nj. A linear expression relating these two was 
demonstrated in several high-energy experiments (II. 5p . 
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< Nb{Nf) >= a + bNf (3.1) 

where the slope, b, is the correlation strength and the intercept, a, is a measure 
of the uncorrelated particles. Both b and a are functions of the energy and atomic 
number of the colliding species. In terms of N^ and N/, a and b can be expressed as, 

<Nb>< Nj > - < NfNb >< Nf> ^ ^ 

- - ^ ^ - — (3.2a) 



< N] > - < Nf >^ 

< NfNb > - < Nf ><Nb> 

< N'j > - < Nf >'^ 



(3.2b) 

^j ^ - <^i\f - - 

Multiplicity fluctuations in high-energy collisions may have three sources: the fluctu- 
ation in the number of particles in each chain (with fixed ends), fluctuations in the 
position of the chain ends (and the subsequent invariant mass of the chain), or the 
fluctuation in the number of chains ^. 

To emphasize long-range correlations, a gap can be introduced between the for- 
ward and backward rapidity windows, eliminating the midrapidity region from con- 
sideration. This has the effect of removing most of the short-range correlations, such 
that predominantly long-range correlations are present. Under the assumption that 
short-range correlations are confined to individual chains, these long-range correla- 
tions are due to the superposition of a fluctuating number of chains. This leads to the 
prediction that the short-range forward-backward correlation strength will decrease 
rapidly with increasing rapidity interval. 

Though the DPM reproduces the results of a great deal of experimental data, 
the assumption of independent chains leads to the realization that even with a large 
density of non-interacting, overlapping strings, they would not provide a means to 
produce a thermalized system. It is predicted that even a small interaction between 
strings will lead to effects that precede the possible formation of a quark-gluon matter 
[6J. 
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3.2.3 Long-Range Correlations in the DPM 

Long-range correlations span a pseudorapidity gap of |?7| > 1.0. The physical 
mechanism that produces long-range correlations is the fluctuation in the number of 
elementary inelastic collisions, controlled by unitarity [12]. The requirement of uni- 
tarity conserves quantum mechanical probability by allowing multiple scattering to 
occur in the model. 

The forward-backward correlation in pp collisions can be expressed as |13] . 

< NfNh > - < Nf >< Nb > = < n > {< NofNoh > - < iVo/ >< iVob >) 



+ 



<n'^ > - <n>'^ 



< Nof >< Nob > (3.3) 



The flrst term in Equation 13.31 corresponds to the correlation between particles 
from a single inelastic collision. These correlations have a short-range in rapidity. 
The second term details the contribution to the long-range correlation. The quantity 
< n^ > — < 72 >^ is the fluctuation in the number of elementary inelastic collisions. 
Introducing a gap about midrapidity (|?7| > 1.0) will substantially reduce the magni- 
tude of the short-range term, which will eventually become negligible with increasing 
1] gap. The average charged particle multiplicity (< A^ >) and average number of 
inelastic collisions {< n >) is given by, 

< A^ >=< n>< No> (3.4a) 

<n>= ^ (3.4b) 

ra=l 

where A'^o is the charged particle multiplicity per inelastic collision and an is the 
probability of n inelastic collisions. In more complicated systems (e.g., collisions of 
heavy nuclei), there are contributions from strings of two varieties, diquark-quark (qq- 
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q) (from the first inelastic collision) and q — q (from subsequent inelastic collisions). 
Then, Equations I3.4al and 13.31 becomes (for large Ar^) 



< A^> ^ 



oo 

71=1 



^ (T„ (2 < iV'?'?-« > +(2n - 2) < iV''-^ > 



.n=l 



(3.5) 



< A^/A/";, > - < A/"/ >< A/";, > = 4 



< n^ > - < n >^ 



< A^"-"? >/< m-'' >b> 

(3.6) 



3.3 Parton String Model 

The introduction of string interaction into the DPM has been accomplished in 
the Parton String Model (PSM) [HHTB] . It is possible that at large string densities 
the assumption of string independence in the DPM may be too strong. If two strings 
overlap, in transverse area, their color fields will overlap, and they may fuse. This 
forms a single string with higher color at the ends. When the fused string hadronizes, 
one will see a reduction in the produced particle multiplicity due to fewer overall 
strings, and an increase in mean p^, a result of energy-momentum conservation. This 
result has been seen experimentally in central high-energy heavy ion collisions. The 
produced particle yield is less than that expected from N^j^ scaling of multiple inde- 
pendent pp collisions. Nun is the number of binary nucleon-nucleon collisions that 
occur when two nuclei collide. These results are also seen in parton saturation phe- 
nomena models, such as the color-glass condensate (CGC) [T7]. The PSM has been 
shown to be in good agreement with pp and heavy ion data at various energies and 
experiments including SPS at CERN and RHIC ^Mi- 
There are two contributions to the multiplicity in heavy ion collisions. The first is 
proportional to Np^rt (the number of nucleons participating in the collision), the sec- 
ond to Nfejji [20]. The experimental results from Au+Au collisions at RHIC are below 
the prediction of Nfej„ scaling from pp to Au+Au. Fusion of strings can account for 
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this multiplicity reduction, which is on the order of ~ 30% in central collisions. The 
PSM also predicts the average transverse momentum {< pr >) enhancement seen in 
Au+Au collisions at RHIC. This is the basic result of overlapping strings conserving 
energy and momentum (fewer particles, but with higher < p^ > per particle). There 
is also good agreement with the inclusive p^ spectra and the pseudorapidity density. 

A predicted consequence of string fusion is the presence of long and short-range 
rapidity correlations. Equation I3.2bl for b can be expressed in terms of the forward- 
backward dispersion squared normalized by the forward (or backward) dispersion 
squared, 

b = 7# (3.7) 

It is predicted that if one studies the long-range correlations by introducing a 
rapidity gap, at large values of Nj, such as those in central heavy ion collisions, b will 
deviate from the expected linear increase [6]. If there is soft (p^ < 2 GeV) string 
fusion, the slope value is expected to decrease even further. Since the number of 
strings increases with energy, atomic number, and centrality, one expects the effects 
of fusion to grow in these cases. This Monte Carlo model (PSM) only considers 
the fusion of two strings close enough in transverse space, and does not fuse hard 
(high-px) strings at all. The model also incorporates a basic rescattering of produced 
particles and decay of resonances. 

3.4 Color Glass Condensate/Glasma 

The Color Glass Condensate (CGC) model is a description of high density gluonic 
matter at small Bjorken x (see Appendix |X]) |21]. In the CGC model, the incoming, 
Lorentz contracted nuclei, are envisioned as two sheets of colliding colored glass. In 
the initial stage of the collision, the color electric and magnetic fields are oriented 
transverse to the beam direction (and each other). After the collision there are 
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additional sources formed that produce longitudinal color electric and magnetic fields. 
This state is referred to as the Glasma, an intermediate phase between the CGC and 
Quark Gluon Plasma (QGP) [22]. This is shown schematically in Figure [37^ |23j . The 
longitudinal fields provide the origin of the long-range correlation and are similar to 
the strings as formulated in the DPM. 

Experimental evidence from the BRAHMS collaboration suggests that the CGC 
is a valid description of the initial state of colliding nuclei in d + Au collisions [24j . 
BRAHMS observes a suppression of the nuclear modification factor Rd+Au (Equation 
II. 4p at high p^ and forward rapidity, in 0-10% most central d + Au events at an energy 
of ^/snn = 200 GeV. This is shown in Figure 13.51 

At midrapidity, for d + Au collisions at an energy of a/stvJv = 200 GeV, there is 
no suppression of Rd+Au, whereas a large suppression is seen in Au+Au at the same 
energy both for mid- and forward rapidities. This indicates that the suppression in 
Au+Au is not due to an initial state effect. The suppression of Rd+Au at forward 
rapidity should be related to the initial conditions of the colliding nuclei, as there is 
no contamination of the final state from a produced quark-gluon medium. 
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(a) Two Lorentz contracted, CGC nuclei incident on one another. The initial 
chromoelectric and magnetic fields are oriented transverse to the direction of 
travel. 





(b) Following the collision, both nuclei recede from one another. Additional 
field sources are induced in each nucleus, creating longitudinal color electric 
and magnetic fields between the nuclei. These fields are similar to the strings 
in the DPM. 



Figure 3.4. Schematic of the CGC/Glasma [23] . 
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Figure 3.5. The nuclear modification factor [R^^au) measured by tlie 
BRAHMS experiment as a function of pseudorapidity (77) in 10% most 
central d + Au collisions at an energy of y^J/wv = 200 GeV. The high-p^ 
suppression at forward pseudorapidity has been suggested as a signal for 
the Color Glass Condensate (CGC) [24j. 
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4. Experimental Analysis 



4.1 Data Analyzed 



Events from the data samples in Table 14.11 were analyzed and the requisite quan- 
tities extracted from the files. This included event-by-event total charged particle 
multiplicities (for centrality determination), charged particle multiplicities in each 
forward and backward measurement interval, and z-vertex position. Every event 
analyzed was from the minimum bias triggered data sample. For heavy nucleus col- 
lisions these were subdivided into centrality bins. Centrality is a characterization of 
heavy ion events into categories based on the measured charged particle multiplicity 
(reference multiplicity). The criteria for determining the reference multiplicity is as 
follows: charged tracks from the primary vertex, fit points in the TPC > 10, within 
the pseudorapidity region —0.5 < rj < 0.5, and with a distance of closest approach 
(dca) of less than 3 cm to the primary vertex. Each centrality bin spans 10% of the 
total multiplicity; from most central to most peripheral: 0-10, 10-20, 20-30, 30-40, 



^s^N (GeV) 


System 


400 


PP 


200 


Au+Au, Cu+Cu, PP 


62.4 


Au-fAu, Cu+Cu, pp 


22.4 


Cu+Cu 



Table 4.1 
A summary of the available data for the various colliding systems and 
energies at RHIC. 
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Figure 4.1. Raw charged particle multiplicity (Nch) in a/s/vw = 200 GeV 
Au+Au showing divisions by centrality. 



40-50, 50-60, 60-70, and 70-80% of the total hadronic cross section. These are the 
centralities used for Au+Au data, while up to 50-60% was considered for Cu+Cu 
collisions. The STAR reference multiplicity distribution at 200 GeV/A in Au+Au 
can be seen in Figure 14.11 

Several quality cuts were implemented during the data processing phase. These 
include some that are also used in the reference multiplicity determination, as men- 
tioned above. The cuts include the number of fit points in the TPC (minimum of 10) 
and dca < 3 cm. 



4.2 Calculating the Forward-Backward Correlation Strength 

The forward-backward (FB) correlation strength b is calculated for all charged 
particles within the STAR TPC acceptance of \r]\ < 1 and with transverse momentum 
Pt > 0.15 GeV/c. The pseudorapidity range is subdivided into forward and backward 
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Centrality 


iVch 


Npart 


Nb^n 


0-10% 


Nch > 431 


325.9 +5.4 -5.3 


940.0 +66.9 -69.5 


10-20% 


312 < Nch < 430 


234.6 +8.3 -9.3 


591.3 +51.9 -59.9 


20-30% 


217 < Nch < 311 


166.7 +9.0 -10.6 


368.6 +41.1 -50.6 


30-40% 


146 < Nch < 216 


115.5 +8.7-11.2 


220.2 +30.0 -38.3 


40-50% 


94 < Nch < 145 


76.6 +8.5 -10.4 


123.4 +22.7 -27.3 


50-60% 


56 < Nch < 93 


47.8 +7.6 -9.5 


63.9 +14.1 -18.9 


60-70% 


30 < Nch < 55 


27.4 +5.5 -7.5 


29.5 +8.2-11.3 


70-80% 


14 < Nch < 29 


14.1 +3.6 -5.0 


12.3 +4.4 -5.2 



Table 4.2 
Centrality definition in terms of the number of charged particles for 
Au+Au at y/SNN = 200 GeV. The number of nucleons participating in the 
collision (Npart) and the number of binary nucleon-nucleon collisions {Nbm) 
estimated from Monte Carlo Glauber simulations |25] are also shown. 
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Figure 4.2. Schematic diagram of the measurement of a forward-backward correla- 
tion. 



measurement intervals of width 0.2 rj. The FB correlation strength b was measured 
in forward and backward intervals that are symmetric about midrapidity [rj = 0). 
The separation between the intervals was measured from the center of each bin and 
included Ar] = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, and 1.8 units in pseudorapidity. 
The bins at At] = 0.2 are contiguous in t] space at ?7 = 0. This is the only bin for 
which there is no physical gap in t] space. 

The coordinate system defined above differs from that of the standard two par- 
ticle pseudorapidity correlation [26]. The FB correlation strength is measured in a 
coordinate system where 77 = is always physically located at midrapidity (Figure 
14. 2p . The collision vertex (z = 0) defines this point. Therefore, all Ar] values are 
measured in an absolute (fixed) coordinate system about i] = 0. In the usual two 
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Figure 4.3. (a) Mean forward charged particle multiplicity (< N ^ >) and 
(b) Mean backward charged particle multiplicity (< N\, >) versus total 
charged particle multiplicity fitted with a linear polynomial up to Nch = 
0-600. (c) < Nf*Nf >) and (d) < Nf*Nh >) versus total charged particle 
multiplicity fitted with a second order polynomial up to Nch = 0-600. 



particle correlation measurement, the relative f] difference between the particles is 
considered. This has particular implications for the definitions of the measurement 
interval and how the centrality determination is made. 

A strong bias exists if the measurement of N^h overlaps with the measurement 
interval in rj space. Therefore, there are three rj regions used to determine reference 
multiplicity: 
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Figure 4.4. The pseudorapidity {rj) distribution of charged particles in 
minimum bias Au+Au data at ^/snn = 200 GeV. 



• The FB strength calculated for At] > 1.0 used the charged particle multiplicity 
(Nch) measured in the STAR TPC r] range from \t]\ < 0.5. 

• The FB strength calculated in the interval Ai] < 1.0 used N^t measured in the 
STAR TPC from 0.5 < \r]\ < 1.0. 

• The FB strength calculated for Ar] = 1.0 was complicated by the fact that the 
measurement window overlapped both of the previously utilized Nch measure- 
ments. Therefore, this particular data point used the charged particle multi- 
plicity measured in the STAR TPC from \ri\ < 0.3 + 0.6 < \r]\ < 0.8. 

This is possible due to the flat t] distribution of the STAR TPC within the mea- 
surement ranges being considered. The t] distribution for charged particles in the 
TPC for Au+Au collisions at ^/sNN = 200 GeV is shown in Figure 14.41 
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For this analysis, the forward (backward) interval is taken as positive (negative) 
7]. In the case of symmetric collisions, these definitions are purely ancillary. 

In order to eliminate the effect of statistical impact parameter (centrality) fluc- 
tuations on the measurement of the FB correlation strength, each relevant quantity 
(the uncorrected mean quantities (A^/),„,„„, (A^6)„„eorr' (^/)Lcorr' and {NfNh)^^^^J 
was obtained on an event-by-event basis as a function of STAR reference multiplicity, 
Nch (Figure lO])- A linear fit to (A^/)„„corr ^^^ {^b)uncorr^ o^ ^ second order poly- 
nomial fit to (Nf) and (NfNf,) , was used to extract the average of these 

\ J ' uncorr \ J ^' uncorr' ^ 

quantities as functions of Nch- Tracking efficiency and acceptance corrections were 
then applied to each event (Section 14.2. ip . Due to statistical limitations, it is not 
possible to apply corrections for every value of Net- One value of the correction, 
calculated for each centrality bin, is applied to (A^/)„„co„. , {Nb) uncorr ' (^/)Lcorr' ^^^ 
(NfNb)^^^^^^ for every event in that centrality bin. Therefore, all events falling within 
a particular centrality have the same efficiency correction. The corrected values of 
(Nf), (Nb), {Nf) , and {NfNf,) were then used to calculate the backward-forward and 
forward- forward dispersions, D^r and D'jr, binned according to the STAR centrality 
definitions in Table 14.21 and normalized by the total number of events in each bin. 

4.2.1 Efficiency Correction 

Because the detectors do not operate perfectly (the tracker may miss particles, 
especially in a central heavy ion events, the tracks may be split or merged, etc.) 
and do not cover a complete volume (e.g., due to gaps between sectors), a correction 
for tracking efficiency and geometrical detector acceptance must be computed. This 
efficiency correction is defined as. 



^^^^ (4.1) 



-' ' sivau 
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where N^ec is the number of reconstructed charged particles and Nsimu is the total 
number of simulated charged particles. When similar cuts are applied to the sim- 
ulated and real tracks, the acceptance effect is implicitly included. The simulated 
particles are produced with the heavy ion event generator, HIJING [27]. These simu- 
lated particles are propagated through a virtual representation of the STAR detector 
constructed in the GEANT simulation framework [28j. The simulated particles inter- 
act with the virtual detector and its constituent materials. Some particles are lost to 
conversions in the material or other mechanisms, while the rest can be found by the 
reconstruction algorithm. The efficiency is a function of collision centrality, particle 
PTi and pseudorapidity. The TPC is least efficient in high density central collisions 
(~ 80%), but improves to over 90% in peripheral events. A sample of the correction 
values is shown in Table H73] for 0-10% ^/s^n = 200 GeV Au+Au collisions. In central 
collisions at this energy the single particle tracking efficiency is ~ 80% while the two 
particle efficiency drops to ^ 65%. 

Forward-Backward Multiplicity Correlation Efficiency 

To accurately determine the forward-backward correlation strength ( b) both single 
particle and two particle tracking efficiency must be determined. From the definition 
of the correlation strength Equation I3.2bt 

< NfNh > - < Nf >< Nh > 
~ < Nj > - < Nf >^ 

it is seen that the single particle efficiency is necessary to correct Nf and A^^. The 

efficiency for finding two particles in the forward pseudorapidity hemisphere (eAr^^Ar^,) 

and one particle in the forward and one in the backward hemispheres {GNf*Nt)y is also 

required. For each of the four quantities, a correction factor is obtained as discussed 

in Section l4.2.1l for each centrality, pseudorapidity interval, z- vertex cut, and number 

of fit points on the track. This correction was applied to Nf, N^, Nf * Nf, and Nf * Nfy 

in every event. 
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Ar^ 


^Nf 


fAffc 


^Ni*Nf, 


^Nf*Nf 


0.2 


0.816 


0.804 


0.651 


0.662 


0.4 


0.825 


0.804 


0.659 


0.678 


0.6 


0.833 


0.809 


0.670 


0.691 


0.8 


0.844 


0.821 


0.689 


0.709 


1.0 


0.854 


0.834 


0.708 


0.726 


1.2 


0.863 


0.844 


0.725 


0.743 


1.4 


0.873 


0.855 


0.744 


0.761 


1.6 


0.873 


0.853 


0.742 


0.761 


1.8 


0.847 


0.821 


0.691 


0.716 



Table 4.3 
Efficiency corrections as a function of Ar] for Au+Au, ^/snn = 200 GeV, 
0-10% centrality, \vz\ < 30 cm, and ffi points > 10. 
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4.2.2 Error Studies 

Statistical and systematic uncertainties have been studied. The error in the mea- 
surement of the forward-backward (FB) correlation strength is predominantly due to 
systematic effects. 

Statistical Errors 

A large sampling of minimum bias data is available for analysis. For central 
Au-|-Au data, the statistical error is less than 1%, while for peripheral it is less than 
3%. This and other contributions from uncertainties in the efficiency estimates are 
taken into account by including an additional 20% in the overall error estimate for 
the FB correlation strength. 

Systematic Errors 

Three main sources of systematic error were considered: the z-vertex position 
with respect to the center of the TPC (3), the fit range (in Net) of the polynomial 
functions to the raw data (3), and the number of TPC fit points on each track (2). 
This provides 3 x 3 x 2 = 18 different combinations that can be used to calculate 
systematic deviations from the mean. The total error on the FB correlation strength 
b is calculated from the errors on Dfr and D'j^ assuming they are uncorrelated, but 
that the error propagated to b is correlated. This is accomplished by calculating the 
standard deviations of Dh and Djr from the 18 separate values. When the final 
calculation of b is made, these errors are treated as correlated since there is overlap 
in the measurements of D^r and D^rr on an event-by-event basis using the same t] 



window in the forward direction. This is shown in Equation 14. 2[ the formula used for 
calculating the correlated errors on b. 
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There is a strong interplay between the three main sources of systematic error as 
a function of Arj. This is predominantly due to the physical limitations of the STAR 
TPC and the tracking algorithm. This has a large effect near the center of the TPC 
and at the outer edges of the acceptance. At small values of At], close to the center 
of the TPC, one track may be split into two if it crosses the central membrane of the 
TPC. At large values of At] there are tracks that may fall outside the acceptance if 
the z-vertex position is substantially shifted from z = 0. The effect on the errors is 
seen as follows (in central Au+Au data): 

1. Changes in the z-vertex position from ± 10 to ± 30 lower the value of b at 
large Ar/ by ~ 15% when the number of fit points on track is 15, but only by 
~ 6% for tracks with 10 fit points. There is a loss of tracks when the z-vertex 
is highly shifted from the center of the TPC and longer tracks (more fit points) 
are required. At small At] for both fit point cuts as a function of Vz, the change 
is less than a few percent. 

2. Changing the fit range affects the 0-10% centrality bin more than the peripheral 
bins. Figure 113] shows that for large values of N^h the fit begins to deviate from 
the data. 

3. Adjusting the number of fit points on the track predominantly affects the FB 
correlation strength at small and large values of Arj. This is related both to the 
position of the z-vertex (for small Arj and due to track loss at the outer edge of 
the measurement {Arj = 1.0)). 
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Average Value 

The final value of the correlation strength for each At] position and centrality is 
the arithmetic mean of the 18 values, determined by varying the cuts as discussed in 
Section 11221 
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5. Experimental Results and Discussion 

The following sections present the results for the forward-backward (FB) correlation 
strength b, calculated as discussed in Chapter HI All plots of D^jr, Djj-, and b are 
corrected for efficiency and acceptance (Section 14.2. ip and shown with the calculated 
statistical and systematic errors (Section I4.2.2p . 

5.1 FB Correlation Strength in Proton-Proton Interactions 

5.1.1 Energy Dependence of the FB Correlation Strength in pp 

The forward-backward (FB) correlation strength in pp collisions has been ana- 
lyzed. The data from pp collisions provides a baseline for the measurement of the FB 
correlation strength that is not contaminated by initial state effects of the colliding 
nuclei or final state effects from the possible production of a quark-gluon matter in 
A+A collisions. Data was obtained at three energies: -x/s^n = 400, 200, and 62.4 
GeV. 

The forward-forward (D'jr) and forward-backward (-D^r) dispersions for ^/SNN = 
62.4 GeV minimum bias pp data are shown in Figure [5TT] as a function of the pseudo- 
rapidity gap. At]. D'jf is shown by the diamond symbols, while D^r are shown by the 
stars. The ratio D^ J D'jf is the FB correlation strength b. Djr is approximately fiat 
as a function of At], while D^f exhibits a decreasing trend. This is refiected in Figure 



15. 2[ the FB correlation strength b for ^/sNN = 62.4 GeV minimum bias pp collisions. 
The FB correlation strength shows only a short-range contribution {At] < 1.0). For 
Ai] > 1.0 the FB correlation strength approaches zero. 
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Figure 5.1. Backward- forward (stars) and forward-forward (diamonds) 
dispersions (D^j and -Djj) as a function of the pseudorapidity gap Ar] in 
minimum bias pp data at ^/sJ^ = 62.4 GeV. 
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Figure 5.2. Forward-backward correlation strength 6 as a function of the 
pseudorapidity gap Arj in minimum bias pp data at yj/viv = 62.4 GeV. 
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Figure 5.3. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (D^j and D}f) as a function of the pseudorapidity gap At] in 
minimum bias ^/sNN = 200 GeV pp data. 



Data from ^/sNN = 200 GeV pp minimum bias collisions for Dfj, Djj, and b are 
shown in Figures [573] and \5M Similar to pp at ^/sJ^ = 62.4, D^^ is flat as a function 
of Arj. At energies as high as y/SNN = 200 GeV, D^j- displays a decreasing trend as a 



function of Ar^. As demonstrated in Figure [531 this directs the behavior of b , which 
shows only a short-range component of the FB correlation strength for Ar/ < 1.0. 

The results for minimum bias pp collisions at ^/sNN = 400 GeV closely resemble 
those at y/sNN = 200 GeV. The FB correlation strength (b) is very similar at both 
energies. At Ar^ = 0, 6 = 0.115 ± 0.002 at ^s^ = 200 GeV and 0.137 ± 0.007 at 
^/sNN = 400 GeV. The FB correlation strength at Arj = is predominantly driven 
by short-range correlations. It is interesting to note that the value of b at large values 
of At] plateaus at the same, small non-zero value for both ^/sNN = 200 and 400 GeV. 
This trend is broken at the lower energy a/J/vat = 62.4 GeV, where b approaches 
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Figure 5.4. Forward-backward correlation strength 6 as a function of the 
pseudorapidity gap Arj in minimum bias ^/sN^ = 200 GeV pp data. 
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Figure 5.5. Forward-backward correlation strength & as a function of the 
pseudorapidity gap /S.rj in minimum bias pp data at ^/snn = 400 GeV. 



zero at large values of Arj. The breakdown of Koba-Nielsen-Olesen (KNO) scaling 
in pp interactions at high energy [291 130] provides a model independent correlation 
between the onset of long-range FB correlations and multiparton interactions in pp 
collisions. The FB correlation strength at At] = is predominantly driven by short- 
range correlations. It is interesting to note that the value of b at large values of 
Arj plateaus at the same, small non-zero value for both y^JJvJv = 200 and 400 GeV. 
This trend is broken at the lower energy ^/snn = 62.4 GeV, where b approaches zero 
at large values of At]. The breakdown of Koba-Nielsen-Olesen (KNO) scaling in pp 
interactions at high energy [29|[30] provides a model independent correlation between 
the onset of long-range FB correlations and multiparton interactions in pp collisions. 
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5.2 Forward-Backward Correlation Strength in Nucleus-Nucleus Colli- 
sions 

5.2.1 Au+Au 

The forward-forward {D'jr) and forward-backward (D^p) dispersions for 0-10% 
^/sNN = 200 GeV Au+Au coUisions are shown in Figure 15.61 as a function of the 
pseudorapidity gap, Arj. D'lr is shown by the diamond symbols, while Dl^ are shown 
by the stars. The ratio D^J D'lr is the FB correlation strength h. Within the total 
error D^jf and D^f is fiat as a function of A?]. This is reflected in Figure [521 the FB 
correlation strength h for 0-10% ^Js^n = 200 GeV Au+Au collisions as a function 
of A?]. The FB correlation strength is also fiat as a function of A?7 within the total 
error, reflecting the behavior of D'lr and Dl^. In the ratio 6, some of the individual 
errors on D']f and D^j are reduced. 

The short-range correlations (SRC) are defined as the correlation between par- 
ticles separated by less than Arj = 1.0. The long-range correlation (LRC) is taken 
as a correlation between particles separated by greater than one unit in At]. Figure 
15.71 shows a strong LRC as a function of Arj. If only short-range correlations (from 
sources such cluster formation, jets, resonance decay, etc.) are present, the expecta- 
tion is a quickly decreasing FB correlation strength as a function of Arj. 

Figures 15.81 and 15.91 show D^jr, Djjr, and b as functions of Ar] for 40-50% ^/snn 
= 200 GeV Au+Au collisions. Unlike in the 0-10% case, the FB correlation strength 
is not flat as a function of Arj, but decreases rapidly as expected in the case where 
predominantly short-range correlations are present, as in pp at an energy of ^/SNN = 
62.4 GeV. From Figure EHl it appears that mid-peripheral Au+Au collisions has the 
expected behavior if only short-range correlations are present and does not exhibit a 
FB long-range correlation. 

Figure 15.81 demonstrates that it is the numerator Dlj that drives the behavior of 
the FB correlation strength. Figure 15.81 shows that D'j^ is approximately flat as a 
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Figure 5.6. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (D^j and Dj-f) as a function of the pseudorapidity gap Ar] in 
central (0-10%) ^/sJ^ = 200 GeV Au+Au data. 
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Figure 5.7. Forward-backward correlation strength & as a function of the 
pseudorapidity gap Ar^ in central (0-10%) i/s/viv = 200 GeV Au+Au data. 
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Figure 5.8. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (-DL and -Dj .) as a function of the pseudorapidity gap Arj in 
mid-peripheral (40-50%) ^/s]^ = 200 GeV Au+Au data. 
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Figure 5.9. Forward-backward correlation strength & as a function of the 
pseudorapidity gap Ar^ in mid-peripheral (40-50%) y/s^^ = 200 GeV 
Au+Au data. 



function of Arj, whereas D^r is maximum at At] = and decreases as a function of 
Ar]. 

The centrality dependence of the FB correlation strength b in Au+Au data at 
^/SlsfN = 200 GeV is shown in Figure 15.101 It is seen that the FB correlation strength 
for central and mid-central (0-10, 10-20, and 20-30%) Au+Au collisions at a/J/vJv 
= 200 GeV exhibits a long-range component that is approximately fiat across Ar]. 
At about the 30-40% centrality bin the long-range component begins to disappear 
and a slight decrease as a function of At] can be seen. By 40-50% it appears that 
all remnants of a long-range correlation have vanished, leaving only the short-range 
component that decreases rapidly with Arj. Figure 15.111 shows the FB correlation 
strength as a function of centrality evaluated at two values of Ai]: Ar] = 0.2 and 1.8. 
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Figure 5.10. Forward-backward correlation strength h for 0-10% (closed 
squares), 10-20% (open squares), 20-30% (closed circles), 30-40% (open 
circles), and 40-50% (closed triangles) most central ^/sNN = 200 GeV 
Au+Au data as a function of the pseudorapidity gap Arj. 



The overall trend shows an increase of the FB correlation strength with centrality, as 
shown in Figures 15.101 and 1 5.111 



Energy Dependence of the Forward-Backward Correlation Strength in 
Au-|-Au 

Additional study of the FB correlation strength was carried out in Au+Au col- 
lisions at a lower energy of ^/snn = 62.4 GeV. As Figures 15.121 and 15.131 show, the 
presence of a long-range correlation in central Au+Au collisions persists at this lower 
energy. For small values of Arj Figure 15.121 shows Dh trending down as a function of 
small Af] (< 1.0), while the opposite is true of D}f- This results in the short-range 
FB correlation strength in Figure 15.131 at At] = having a value close to that shown 
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Figure 5.11. Centrality (Npart) dependence of the forward-backward cor- 
relation strength b evaluated at two At] values: At] = 0.2 (closed squares) 
and 1.8 (open circles). The long-range correlation strength (At^ = 1.8) 
reaches the value of the short-range correlation strength {Arj = 0.2) for 
central collisions {Npart ~ 350). 
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Figure 5.12. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (D^j and D}f) as a function of the pseudorapidity gap At] in 
central (0-10%) ^/s]^ = 62.4 GeV Au+Au data. 



in Au+Au at an energy of y/s^N = 200 GeV (Figure 15. 7p . This is also a result of 
the tracking limitations as discussed in Section 14.2.21 The long-range FB correlation 
strength is lower by approximately 30% at ^/sNN = 62.4 GeV (Figure IB. 13p compared 
to ^/sNN = 200 GeV (Figure ISTTl) . Though the ratio of the energies is approximately 
a factor of three, the ratio of the multiplicities differ by about a factor of 1.5. The 
normalization of b by D'jf contributes to this, since D'jr goes approximately with the 
multiplicity. Therefore, it is the quantity D^r that carries the dynamical information. 
This approximate scaling with multiplicity does not hold when comparing different 
colhding systems (e.g., Au+Au and Cu+Cu). 

For comparison purposes, the 40-50% centrality bin for ^/snn = 62.4 GeV Au+Au 
is shown in Figure l5.14[ As seen in the ^/snn = 200 GeV Au+Au result, the mid- 
peripheral data does not exhibit a long-range component. 
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Figure 5.13. Forward-backward correlation strength 6 as a function of 
the pseudorapidity gap At] in central (0-10%) ^/sNN = 62.4 GeV Au+Au 
data. 
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Figure 5.14. Forward-backward correlation strength 6 as a function of the 
pseudorapidity gap Ar^ in mid-peripheral (40-50%) ^/sNN = 62.4 GeV 
Au-|-Au data. 
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Figure 5.15. Forward-backward correlation strength b for 0-10% (closed 
squares), 10-20% (open squares), 20-30% (closed circles), 30-40% (open 
circles), and 40-50% (closed triangles) most central ^/sNN = 62.4 GeV 
Au+Au data as a function of the pseudorapidity gap Arj. 



The centrality dependence of the FB correlation strength at ^/sNN = 62.4 GeV 
(Figure I5.15P demonstrates the progressive decrease in the long-range component as 
a function of centrality. Additionally, the FB correlation strength in central collisions 
exhibits a dependence on the collision energy. The long-range FB correlation strength 
is smaller at lower energies. 

5.2.2 Cu+Cu 



To supplement the studies of Au+Au collisions at RHIC, collisions of the lighter 
nuclei Cu+Cu were examined at three energies: i/s/vjv = 200, 62.4, and 22.4 GeV. 
The STAR centrality definitions for ^/snn = 200 Cu+Cu data are shown in Table 

KB 
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Centrality 


Nch 


Npart 


-^ ^ bin 


0-10% 


Nch > 140 


98.4 +1.0-1.0 


185.7+5.9-5.2 


10-20% 


103 < N^h < 139 


74.8 +2.5 -2.2 


126.7+6.7-7.1 


20-30% 


74 < Nch < 102 


54.4 +2.8 -2.5 


81.5 +6.0-6.1 


30-40% 


53 < Nch < 73 


38.5 +2.5-3.1 


51.0 +4.8 -5.4 


40-50% 


37 < Nch < 52 


26.3 +2.4-3.3 


30.6 +3.9 -4.6 


50-60% 


25 < Nch < 36 


17.6 +2.6-3.1 


18.2 +3.5-3.4 



Table 5.1 
Centrality definition in terms of the number of charged particles for 
Cu+Cu at ^/sNN = 200 GeV (P06ib production). The number of nu- 
cleons participating in the collision (Npart) and the number of binary 
nucleon-nucleon collisions (Nbin) estimated from Monte Carlo Glauber 
simulations EH] are also shown. 
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Figure 5.16. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (D^j and D}f) as a function of the pseudorapidity gap At] in 
central (0-10%) Cu+Cu data at ^s^ = 200 GeV. 



The study of the FB correlation strength b for Cu+Cu was accomplished in the 
same manner as that for Au+Au. Figure 15.161 shows D^r and D^r for 0-10% most 
central ^/snn = 200 GeV Cu+Cu collisions. Both quantities are fiat as a function 
of Arj, similar to central Au+Au collisions at the same energy (Figure 15^ . Though 
the qualitative trend is the same, the values of DL and D'jr are both much smaller 
in central Cu+Cu than in Au+Au. 

Figure 15.171 is the FB correlation strength b for the 0-10% most central ^/snn 
= 200 GeV Cu+Cu collisions. As with Dh and Dfj, the qualitative trend of b is 
the same for Cu+Cu and Au+Au at the same energy and centrality. Both exhibit 
a strong, long-range component to the FB correlation strength. The plateau of the 
FB correlation strength is ~ 15% lower in Cu+Cu than Au+Au at the same energy 
(Figure [5?7|) . while the multiplicity is ~ 70% lower in Cu+Cu (< N^h >Au+Au^ 116 
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Figure 5.17. Forward-backward correlation strength 6 as a function of the 
pseudorapidity gap Arj in central (0-10%) Cu+Cu data (circles) at ^/sNN 
= 200 GeV. For comparison, the same quantity is plotted for central (0- 
10%) Au-|-Au data at the same energy (squares). 



and < Nch >cu+Cu^ 34). The average multiplicity in 0-10% Cu+Cu data at a/J^tv = 
200 GeV is most closely matched by the 30-40% centrality bin of Au+Au. Comparison 
to the 30-40% Au+Au centrality from Figure 15.101 shows a large difference between 
the FB correlation strength of the two systems at the same average multiplicity. The 
comparison of the FB correlation strength in 0-10% Cu+Cu and 30-40% Au+Au is 
shown in Figure IS. 181 

The FB correlation strength in mid-peripheral ^/snn = 200 GeV Cu+Cu is also 
very similar to Au+Au at the same energy. Figure 15.201 shows the Cu+Cu results 
for the 40-50% centrality bin. The FB correlation strength in mid-peripheral Cu+Cu 
and Au+Au (Figure 15.91) both show the presence of only short-range correlations. 
Both systems show similar qualitative and quantitative evolution as a function of 
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Figure 5.18. Comparison of the forward-backward correlation strength b 
as a function of the pseudorapidity gap Ar] in central (0-10%) Cu-|-Cu 
(squares) and mid-peripheral (30-40%) Au-|-Au (Au-|-Au) data at an en- 
ergy of ^/SNN = 200 GeV. The two centralities have approximately the 
same average multiplicity, but exhibit a different evolution as a function 
of Arj. 
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Figure 5.19. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (D^j and D}f) as a function of the pseudorapidity gap At] in 
mid-peripheral (40-50%) Cu+Cu data at ^/s^ = 200 GeV. 



Arj. There is also a close similarity to the FB correlation strength in pp at the same 
energy. 

The centrality dependence of the FB correlation strength from y/SNN = 200 GeV, 
0-10% to 40-50% Cu+Cu collisions is shown in Figure 15.221 Though the total average 
multiplicity differs dramatically from Au+Au to Cu+Cu, the qualitative trend of the 
FB strength as a function of centrality is similar. The FB correlation strength for 
40-50% Au+Au, Cu+Cu, and minimum bias pp data at an energy of ^/sNN = 200 
GeV are shown together in Figure 15.211 The peripheral heavy ion data exhibits a 
short-range correlation as a function of At], with no long-range component. This is 
in good agreement with pp, which is dominated by short-range correlations. 
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Figure 5.20. Forward-backward correlation strength & as a function of the 
pseudorapidity gap Arj in mid-peripheral (40-50%) Cu-|-Cu data at ^/sN^ 
= 200 GeV. 
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Figure 5.21. Forward-backward correlation strength 6 as a function of 
the pseudorapidity gap Ar] in mid-peripheral (40-50%) Au+Au (squares), 
Cu-|-Cu (triangles), and minimum bias pp (circles) data at ^/snn = 200 
GeV. There is a general lack of long-range correlations in peripheral heavy 
ion data that is consistent with basic proton-proton interactions. 
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Figure 5.22. Forward-backward correlation strength b for 0-10% (closed 
squares), 10-20% (open squares), 20-30% (closed circles), 30-40% (open 
circles), and 40-50% (closed triangles) most central Cu-|-Cu data at ^/sNN 
= 200 GeV as a function of the pseudorapidity gap At]. 
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Figure 5.23. Forward-backward correlation strength 6 as a function of the 
pseudorapidity gap Ar^ in central (0-10%) Cu+Cu data at ^/sJ^ = 62.4 
GeV. 



Energy Dependence of the Forward-Backward Correlation Strength in 
Cu+Cu 

The FB strength was also studied for Cu+Cu collisions at a/J/vtv = 62.4 GeV. The 
mean multiplicity in central (0-10%) Cu+Cu at ^/snn = 62.4 GeV is approximately 
equivalent to that from 10-20% Cu+Cu at ^/s]^ = 200 GeV and 40-50% y/s]^ = 
200 GeV Au+Au. As Figure [5:231 shows for central ^5^ = 62.4 GeV Cu+Cu, the 
FB correlation strength b is about 20% lower than 10-20% y/s]^ = 200 GeV Cu+Cu 
(Figure I5.22p and over four times larger than that from 40-50% a/saT/v = 200 GeV 
Au+Au (Figure ISlOj) . Quantitatively, the long-range component exhibits a similarity 
to both 20-30% Cu+Cu and Au+Au at ^/JJ^ = 200 GeV. 
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Figure 5.24. Forward-backward correlation strength 6 as a function of the 
pseudorapidity gap Ar/ in mid-peripheral (40-50%) Cu+Cu data at a/s/vJv 
= 62.4 GeV. 



The results for mid-peripheral (40-50% centrality) ^/sNN = 62.4 GeV Cu+Cu are 
shown in Figure 15.241 As seen in the same centrality for both Cu+Cu and Au+Au at 
■s/snn = 200 and 62.4 GeV and y/s^N = 200 GeV MB pp, the FB correlation strength 
appears to possess only a short-range component that decreases with increasing Arj. 
The behavior as a function of centrality is shown in Figure 15.251 

Additionally ,the FB correlation strength in Cu+Cu collisions was studied at the 
lowest available energy, y/s^iy = 22.4 GeV. Figure 15.261 shows the result for the 
most central (0-10%) data. The FB correlation strength b is both qualitatively and 
quantitatively similar in central ^/sNN = 22.4 and 62.4 GeV Cu+Cu collisions (Figure 
I5.23p . Though the mean multiplicities are ~ 30% different, the similarity may indicate 
that the particle production mechanism for Cu+Cu at ^/SNN = 22.4 GeV is the same 
as that at ^/sNN = 62.4 GeV. 
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Figure 5.25. Forward-backward correlation strength b for 0-10% (closed 
squares), 10-20% (open squares), 20-30% (closed circles), 30-40% (open 
circles), and 40-50% (closed triangles) most central Cu+Cu data at ^/s^n 
= 62.4 GeV as a function of the pseudorapidity gap Ar]. 
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Figure 5.26. Forward-backward correlation strength & as a function of the 
pseudorapidity gap Ar^ in central (0-10%) Cu-|-Cu data at ^/sNN = 22.4 
GeV. 



76 



cvj K 



CM n 

o 



7 

6 

5 

4 

3 

2b 

1 



"^0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Ari 



Figure 5.27. Backward-forward (stars) and forward-forward (diamonds) 
dispersions (D^j and D}f) as a function of the pseudorapidity gap At] in 
mid-peripheral (40-50%) Cu+Cu data at y/sj^ = 22.4 GeV. 
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Mid-peripheral (40-50%) y/sj^ = 22 A Cu+Cu results (Figure E2S]) is also in close 
qualitative and quantitative agreement with the corresponding measurement at ^/sNN 
= 62 A GeV (Figure [5. 24p . This also holds for the remaining centralities at ^/SNN = 
22 A (Figure [5. 29p and ^/sNN = 62.4 (Figure l5.25p . These are the first indications of 
the energy independence of the FB correlation strength in Cu+Cu below a collision 
energy of y/s^N = 62.4. 



5.3 Discussion 

The study of long-range FB correlations in Au+Au, Cu+Cu, and pp reveals de- 
tails of particle production for the bulk of particles in these collisions. Whether 
quark-gluon matter is produced or not, the vast majority of particles have transverse 



77 



0.7- 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

-0.1 



22.4 GeV Cu+Cu, 40-50% 



I I I I I I I I I I I I I I I I I I I I I I I I I I . I I . I I I 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Ar| 



Figure 5.28. Forward-backward correlation strength & as a function of the 
pseudorapidity gap Ar^ in mid-peripheral (40-50%) Cu-|-Cu data at ys/vJv 
= 22.4 GeV. 
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Figure 5.29. Forward-backward correlation strength b for 0-10% (closed 
squares), 10-20% (open squares), 20-30% (closed circles), 30-40% (open 
circles), and 40-50% (closed triangles) most central Cu+Cu data at ^/s^n 
= 22.4 GeV as a function of the pseudorapidity gap Ar]. 
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momentum, pt < 1-0 GeV. In this analysis, all charged particles with pt > 0.15 
GeV are considered. This probes the behavior of the bulk material formed in these 
experiments. 

The Dual Parton Model (DPM) ascribes the long-range correlation to fluctuations 
in the number of elementary, inelastic collisions, directed by unitarity (Section I3.2.3P 
|12] . In the DPM, the single particle inclusive spectrum in elementary pp collisions is 
expressed as [31j . 



dNPP 



(y) = T.-T.^n {K'^'^iy) + N^r''{y) = i^^n - 2)Nr'<y)) (5.1a) 



dy 

^N^-'^iy) + Nr'\y) + {2k - 2)Nr''iy) (5.1b) 

where k is the average number of inelastic collisions (Equation 13. 4b|) . Each term 
contributes two strings per inelastic collision. Two strings are of the type diquark- 
quark, which come from the valence quarks of the proton. Any additional strings 
manifest between q — q from the quark sea. Recalling Equation 13.31 for the FB corre- 
lation in pp, 

< NfNb > - < Nf >< Nb > = < n > {< NofNob > - < A^o/ >< Not >) 



+ 



< n^ > - < n >^ 



< Nof >< Nob > (5.2) 



with the first term in Equation 13.31 the FB correlation of particles produced in 
a single elementary inelastic collision, while the second term describes the fluctua- 
tion in the number of elementary inelastic collisions, identifled with the long-range 
correlation. 

The generalized version of Equation I5.1al for nucleus-nucleus collisions at flxed 
impact parameter (b) is written as 
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^—(6) = n^ib) [Nll-'"{y) + iV;^,^!^) + {2k - 2)iV;(,7=]+(n(6) - n^{h)) 2kNl:^'^{y) 

(5.3) 

In Equation 15. 3[ ua is the average number of participating (wounded) nucleons in 
the colhsion, n is the average number of binary nucleon-nucleon colhsions, and /i(6) 
is the average number of inelastic colhsions by a particular nucleon at fixed impact 
parameter. The first term in Equation 15.31 describes the nucleon-nucleon interactions 
outlined in Equation I5.1bt scaled by the average number of participating nucleons, 
TiA- In this case, the average number of inelastic collisions is fc ■ n, making the total 
number of strings 2k ■ n. The second term is required at high energies to account for 
the increasing influence oi q — q sea strings [32] . 

It was shown that the FB correlation strength for pp and mid-peripheral (40-50%) 
heavy ion data is dominated by short-range correlations. There is good qualitative 
and quantitative agreement for all energies and colliding systems when considering 
these mid-peripheral results. The lack of long-range FB correlations for peripheral 
heavy ion events indicates that there is a lack of multiple parton scattering. There 
is a strong suggestion that the correlation between particles in mid- and peripheral 
heavy ion collisions is very similar to that in pp collisions, where multiple parton 
interactions are expected to be small. 

It is possible that a small, residual long-range correlation is present in high energy 
pp collisions, perhaps related to the breakdown of KNO scaling above energies of 
■,/snn = 100 GeV [2S1ES]- It has been suggested that this breakdown in KNO scaling 
is related to the increasing contribution of multiple parton interactions in pp and 
pp at higher energies [30]. This provides a model independent correlation between 
long-range FB correlations and multiparton interactions in pp collisions. 

In contrast, central and semi-central ^/sNN = 200 GeV Au+Au and Cu+Cu col- 
lisions exhibit a large FB correlation strength for Ar/ > 1.0. There is a substantial 
difference between the long-range FB correlation strength in central nucleus-nucleus 
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Figure 5.30. Comparison of the measured forward-backward correlation 
strength & as a function of the pseudorapidity gap Arj in central (0-10%) 
Au+Au data at ^/sNN = 200 GeV (closed circles) to the prediction of 
the Parton String Model (PSM) (open circles). The PSM reproduces the 
qualitative trend of the long-range correlation, but under predicts the 
magnitude. For this simulation, the string fusion option is not enabled. 



collisions compared to the pp-like behavior of mid-peripheral and peripheral data. 
The experimental measurement of substantial long-range correlations is indicative of 
multiple partonic interactions in central and semi-central nucleus-nucleus collisions. 



Figure IS. 301 demonstrates the comparison to the Parton String Model (PSM, Sec- 
tion [S3]) for central 0-10% Au+Au at an energy of ^/sj^ = 200 GeV. The PSM 
reproduces the qualitative trend of the FB correlation strength as a function of Arj. 
This provides a model dependent interpretation of the long-range FB correlation as 
being due to multiple partonic interactions. 
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A possible qualitative explanation for the centrality dependence of the FB correla- 
tion strength as shown in Figure \5.10\ relates to the idea of the formation of a central, 
partonic core surrounded by a hadronic corona. This has been discussed in terms of 
hydrodynamic properties and finite formation time of the strongly interacting quark 
gluon plasma (sQGP) |Mll35], as well as a model to specifically quantify individual 
contributions from the core and corona ^U\. The primary difficulty is to find an 
experimental signal that can distinguish the respective contributions of the core and 
corona. The measurements of various observables integrate over all particles produced 
in the collision, whether they are produced from a single source (i.e., partonic system) 
or typical pp-like (i.e., hadronic) interactions. If the manifestation of a large, long- 
range correlation is due to partonic effects (as described in the DPM), the partonic 
core formed in heavy ion collisions would be the primary source of this correlation. 
Since the hadronic corona would produce particles similar to pp interactions (with no 
formation of partonic matter), the contribution of the corona would consist of only 
short-range correlations. It is expected from geometrical arguments that this partonic 
core would have a larger volume in central heavy ion collisions, compared to more 
peripheral interactions. Therefore, the increasing infiuence of a partonic core would 
lead to a larger long-range correlation. The evolution with centrality also agrees with 
predictions from the DPM pTHH] . 
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6. Summary 

6.1 Summary 

The energy and system size dependence of long-range multiplicity correlations has 
been studied using the forward-backward (FB) correlation strength b. The FB corre- 
lation strength was calculated in both nucleus-nucleus and hadron-hadron collisions 
for three systems: Au+Au, Cu+Cu, and pp. The energies considered were: ^/SNN 
= 200 and 62.4 GeV for Au+Au, ^s^ = 200, 62.4, and 22.4 GeV for Cu+Cu, and 
^ysNN = 400, 200 and 62.4 GeV for pp. In the most central nucleus- nucleus collisions 
at ^ysNN = 200 GeV, the magnitude of the FB correlation strength is approximately 
flat across a wide range in Arj. This includes a strong, long-range correlation for 
At] > 1.0. The long-range multiplicity correlations in central nucleus-nucleus colli- 
sions appear to depend less on the size of the colliding system (Au+Au or Cu+Cu) 
than on the energy of the incident nuclei. The data demonstrates that the long-range 
FB correlation strength in central heavy ion collisions is predominantly driven by 
the collision energy. This is the first indication of dense matter, which also exhibits 
partonic characteristics, in central heavy ion collisions. 

The long-range correlation decreases with decreasing centrality. By ^ 40-50% 
most central Au+Au (and Cu+Cu) collisions, the FB correlation strength is con- 
sistent with only short-range correlations, as also observed in pp interactions at the 
same energy. The evolution as a function of At] in semi-peripheral data shows a 
similar behavior across all energies and colliding systems. The energy dependence of 
the FB correlation strength in pp shows that for ^/SNN = 200 and 400 GeV, the FB 
correlation strength plateaus at a similar (small) value for Arj > 1.0, while y^JJvJv = 
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62.4 GeV data goes smoothly toward b = 0. 

These long-range correlations can be ascribed to multiple elementary inelastic col- 
lisions, which are predicted in the Dual Parton Model and the Color Glass Conden- 
sate/Glasma phenomenology. Additionally, the association of long-range correlations 
with the breakdown of KNO scaling provides a model independent interpretation in 
terms of multiple partonic collisions. The centrality dependence of the FB correlation 
strength demonstrates the increasing dominance of a partonic core in the resulting 
particle production. This indicates that substantial amounts of dense partonic matter 
are formed in central Au+Au (and possibly Cu+Cu) collisions at ^/sNN = 200 GeV. 

6.2 Future Study 

Additional studies of the FB correlation strength remains to be done. Some of 
the items requiring further study are: 

• Quantifying the number of particles produced from the core versus the corona. 

• The dependence of the FB correlation strength on transverse momentum (pr)- 

• The particle species dependence of the FB correlation strength, in terms of the 
difference between mesons and baryons. 
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A. Appendix 

The kinematic variables used in relativistic heavy ion (and other high-energy) physics 
are defined in such a way that they have simple properties under Lorentz transfor- 
mation. All notation makes use of natural units, c = h = 1. The rapidity variable, y, 
is defined as follows, 

y = - ('-^^ (A.i) 

This dimensionless quantity is defined in terms of po? the energy of the particle, 
and Pz) the longitudinal momentum of the particle. The longitudinal direction is 
the direction defined by the direction of the beam. An additive constant relates the 
rapidity of particles in different frames of reference. The relationship between p^, p^, 
and y is, 

Po = rriT cosh y (A. 2a) 

Pz = rriT sinhy (A. 2b) 

where the transverse mass, rri^ = m^ +p\,. The rest mass of the particle is m, and 
Pt is the transverse momentum, the component of the particle momentum orthogonal 
to the beam direction. 

To measure rapidity, an experiment must measure both energy and longitudinal 
momentum. A more tractable experimental measure is often the angle of the particle 
with respect to the beam axis. In that case, the pseudorapidity variable is useful. 



In 



tan - 
\2. 



(A.3) 



where 9 is the angle between particle momentum and the beam axis. Pseudora- 
pidity and rapidity become approximately equal for large particle momenta. 

In the collision of projectile A with target B, A + i? — t- A + i?, the collision is 
referred to as elastic. If the target dissociates, A + B ^ A + X, the reaction is called 
inelastic. The momentum transfer is defined as, 

q = Ak = k-k' (A.4) 

where k and k' are the initial and final 4-momenta of the projectile, respectively. 
In inclusive deep inelastic scattering experiments with electrons and protons, the 
scattering cross section was found to be independent of q^ and scaled with the ratio 



X = ^ ^ = , , A. 5 

2q-p 2Miy ^ ' 

This was unlike the elastic cross-section, which is highly dependent on q^. The 
Bjorken x is the fraction of the momentum carried by a parton in the infinite mo- 
mentum frame and v is the energy transfer E — E' . 
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B. Appendix 
B.l FTPC 
B.1.1 Overview 

The Forward Time Projection Chamber (FTPC) subsystem at the STAR experi- 
ment consists of two cyhndrical drift chambers 75 centimeters in diameter and 120 cen- 
timeters in length, mounted at ± 2.34 meters from the center of the Time Projection 
Chamber (TPC). These detectors provide full azimuthal (27r) charged particle track- 
ing and momentum determination over a pseudorapidity [r]) range of 2.5 < \r]\ < 4.0. 
Due to the particle production rates in heavy ion collisions such as Au+Au, as well 
as space constraints within the TPC, the FTPC makes use of a radial drift field to 
maximize two-track resolution. This is unlike the main TPC, which utilizes a longi- 
tudinal electric field to drift electrons to the end caps. The structure of the FTPC 
is shown in Figure IB. II [T]. The readout chambers are curved structures located in 
five rings around the outer surface. These rings consist of two padrows each and are 
divided into six azimuthally equivalent readout chambers. The field cage consists of 
the inner, metallized plastic cylinder which functions as the HV-electrode, the outer 
cylinder wall at ground potential, and a concentric ringed, aluminum field cage that 
closes the field at either end. 

Because the FTPC utilizes a radial field, the drift electrons do not travel parallel 
to the STAR magnetic field as they do in the main TPC. Instead, they drift or- 
thogonal to the magnetic field and experience a Lorentz force proportional to ExB . 
This effect must be taken into account when reconstructing hit positions. At the 
inner radius, where the cluster density and the drift distance are maximal, the Exn 
effect leads to a broadening of the cluster distribution, improving the two-track reso- 
lution capability. The drift electrons are produced by ionization due to the passage of 
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Figure B.l. Schematic of a Forward Time Projection Chamber (FTPC). 
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charged particles through the FTPC gas volume. The gas used is a 50%-50% mixture 
of Ar-C02, selected for its drift characteristics and chemical stability. The FTPC 
readout chambers are curved in order to maintain as close to a radial field as possi- 
ble. Each FTPC has 9600 pads with 256 timebins/pad. Due to space constraints, the 
number of padrows in the FTPC is 10, which is also the maximum number of hits 
on track. Therefore, the FTPC momentum resolution is limited to ~ 20% across the 
range 0.1-2.0 GeV/c. Energy loss calculations, dE/dx, are also not possible with the 
FTPC due to the limited number of fit points. Therefore, the FTPC does not have 
particle ID capability. 



B.1.2 FTPC Calibration 

The information presented in this section consists of several steps that are neces- 
sary to maintain the FTPC subsystem. In the past several years, the various software 
and hardware tasks to calibrate the FTPC for data taking purposes have been learned 
and implemented. This has required extensive time at Brookhaven National Labora- 
tory in support of FTPC operations. All these tasks needed to be distilled down to 
the point where they could be carried out by one student. Previously, the support 
for the FTPC consisted of a larger group from another STAR institution. A recently 
completed, extensive write-up describing FTPC operational parameters and calibra- 
tion procedures will help ease the training of additional people in FTPC procedures. 
This work is considered service work for the STAR collaboration. 

B.1.3 FTPC Software and Operational Tasks 

The FTPC software is an extensive set of C++ programs integrated into the 
STAR computing environment. The major programs can be broken down into four 
general categories: 
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1. Cluster finding; 

2. Particle tracking; 

3. Laser calibration; 

4. Diagnostic/operations. 

The FTPC must be calibrated for every run. Additionally, if the collider energy 
or particle species is changed, several parameters may need to be adjusted. There 
are several required procedures to both calibrate the detector and maintain detector 
effectiveness during and after data taking, including: 

1. Laser analysis to verify gas composition; 

2. Inner cathode correction of the internal detector geometry; 

3. Transverse (x, y) FTPC vertex offset with respect to the TPC, which takes into 
account rotational corrections about FTPC mounting points; 

4. Longitudinal (z) FTPC vertex offset with respect to the TPC, which indicates 
changes in electronics timing (to); 

5. Analysis procedure for determining effective inner radius for cluster reconstruc- 
tion (cluster radial step position), which provides a check on measured temper- 
atures and to; 

6. Gain tables to mask dead/noisy front-end electronics (FEEs); 

7. Embedding to verify the tracking efficiency of the FTPCs; 

8. Working knowledge of the FTPC control software, including interactive graphical 
user interfaces (GUIs) that are used to directly monitor and control the real-time 
operation of the detector; 

9. Continual monitoring of quality assurance (QA) plots to diagnose potential prob- 
lems during several months of data taking; 

10. Troubleshooting of problems that arise both during the run and in fully recon- 
structed data. 
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Figure B.2. Laser configuration in one FTPC. 



B.1.4 FTPC Laser System 



The gas composition of tlie FTPC determines tlie drift characteristics of the de- 
tector. Therefore, as part of the FTPC cahbration process, it must be determined 
whether the gas composition has changed over the course of the data run. The FTPC 
laser system is integral for checking the gas composition and ExB corrections. The 
FTPC is sensitive enough to see changes in the gas composition on the order of 1%. 
There are a maximum of 15 laser tracks per FTPC, 5 each in three laser sectors. 
In each sector are 3 straight tracks that run parallel to the beam pipe at a partic- 
ular radius, and 2 inclined (diagonal) tracks. The measured radial positions of the 
3 straight tracks are 11.91 cm, 19.55 cm, and 28.56 cm from the beam line. Laser 
runs are taken every few days during the run, or as needed. The events are processed 
through a modified version of the FTPC reconstruction chain that is optimized for 
laser runs. Once reconstructed, the data is processed once again to plot the positions 
and residuals for the laser tracks, after corrections have been made for magnetic field, 
to, and gas composition. The laser configuration and an example of the data from a 



representative laser run are shown in Figures [6.21 B.3(a) and B.3(b) 



The straight beams are used to check the gas composition and to- Because the 
inner beam has the longest drift distance, it will be more affected by changes to the 
gas composition. Conversely, the outer beam will not be as affected by the gas com- 



94 



I ^3H^tu5(s^^a^ght^ laser Iracks ] 



Mean ».W 
RMS B.065 



1 Numt»rhtt2on1a2.rt«ck | 




nhlts 




Eiiliiis 434 


, 

290 


r 








Mean HJ2 


18D 


r 








RM3 iai? 


ISO 










lints i-^" 


r 












120 














WO 














80 


r 












ea 


r 












40 


H 












ZO 


1- 














P -...(.-. it:-. . H .... 1 . . 


. . 1 . . 


. . 1 . . 


. . 1 . 











5 


10 


15 


20 


25 


30 
















r {cm) 




1 Tadlusfst 


olghl) 


laser tracks 


J 












3D 






. . - 


- 




- - 


Mown J 


2S 






position 
(cm) ^s 


_ 










' * 




CO 


10 

s 


- 
















M 





1S!1 


' '<»' ' JIB 


'zzo' 


%' ' 


3I0' ' 





5 6 7 9 3 to 11 

f^ liils on 1racl< 



I phi iaaer tracKs 1 



500 
400 - 
90 
3Q(»- 



phi 



EnliiBB 3028 
Mean -74.6 
RMS Dl6309 



z (cm) 



^« 50 IW 1» 200 250 

Piii (tieyrees) 



(a) Clockwise from upper left: Laser ADC vs. radial position in the FTPC; # hits 
on track; Phi position; Laser radial position vs. z-position. 



[ X Residual | 



1400 


r 








p 










1200 


- 




















woo 


- 




















mo 


1- 




















eoD 


- 




















400 


r 




















200 


- 




















« 






















j'-o'.t 


-03 


-OJ 


Jl 


■ 


0.1 


oi 


OJ 


0.H 


as 



I r Residual | 





115 -04 -OJ -Oi -0.1 0.1 Oi DJ 0,4 OJ 



I Phi Residual | 




0.5 -04 -0.3 -02 -0.1 0.1 0-Z 03 04 0^ 



(b) Clockwise from upper left: x, y, r, and phi laser track residuals 



Figure B.3. Results from a laser run showing all three radial beams in one 
sector. 
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position, but will be influenced more by changes to to- Ideally, the positions of the 
straight tracks will agree (or be close to) the theoretical positions with no change in 
either to or gas composition. If that is not the case, iterating through various values 
of to and gas composition is required. The residuals of the inclined tracks are used 
to check the Exn corrections. 

The high radiation operational environment in heavy ion collisions takes a physical 
toll on detectors. Occasionally, a large amount of stray particle flux from degrading 
beam conditions or the heavy ion collisions themselves can create currents in electron- 
ics beyond the level deemed safe. The affected anodes "trip", or discharge all their 
energy rapidly to hopefully prevent damage to current sensitive electronics. This is 
not always successful, and damage to Front End Electronics (FEEs), or even entire 
Readout Boards (RDOs), is possible. Dead and noisy electronics reduce FTPC track- 
ing efficiency. Dead areas can be corrected with efficiency results from embedding, 
but noisy electronics are also useless when reconstructing data. The solution is to 
mask out these noisy electronics after data taking, but before data reconstruction. 
This is accomplished using a gain table. The gain table multiplies each pad by a 
calculated gain factor. If a pad exceeds a user defined noise cut (ADC count), the 
gain factor is set to and the noisy pad is now a "dead" region. 

There are additional steps that must be taken to calibrate the FTPCs to produce 
physics data. The FTPCs use tracks to independently reconstruct the primary vertex. 
These FTPC vertices will not be identical to the TPC vertex due to several factors. 
This includes changes in to (which affects the z-component of the vertex position), the 
long lever arm of FTPC tracks to the primary vertex (small errors propagated over 
large distances), and a slight physical shift (or rotation) about the FTPC mounting 
points. This slight shift is exacerbated by the long lever arm when projecting tracks 
to the primary vertex, resulting in an offset of several millimeters in the transverse (x, 
y) plane. This offset needs to be corrected to ensure that the FTPC and TPC vertices 
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match as closely as possible for proper data reconstruction. In 2004 it was discovered 
that the reconstructed FTPC vertices were different for various STAR magnetic field 
settings. There appeared to be an almost independent, orthogonal shift in the ver- 
tex position for each FTPC (one moved mostly in the x-direction, the other in the 
y-direction). When the magnetic field polarity was returned to its initial setting, the 
reconstructed FTPC vertices almost, but not quite, returned to their initial values. 
This could possibly be due to a small movement of the entire STAR detector and 
magnet coils due to a change in the force vector after a magnet polarity flip. 

The usable, inner volume of the FTPC begins at approximately 7.80 cm from the 
beam line. Therefore, good, reconstructed clusters should not be found at distances 
smaller than 7.80 cm. There are a small amount of bad clusters, either from electron- 
ics noise, beam background, or other sources, that are reconstructed below 7.80 cm. 
At ^ 7.80 cm, there is a rapid rise in the cluster count, seen in Figure [K4l 

As previously mentioned, the FTPC utilizes a cylindrical cathode located at the 
inner radius to produce the drift field. This would ideally be located in the direct 
center of the FTPC and provide perfect cylindrical symmetry for the produced electric 
field. However, due to slight machining errors, the cathode is not perfectly centered 
in the FTPC. This shift is on the order of 0.25 mm, but corresponds to a factor 
of 10 increase at the outer radius, or 2 mm. Additionally, the detector is sensitive 
enough that the effect of gravity on the inner cathode, causing a slight warping in the 
vertical direction, is also noticeable. The inner cathode offset has two, sympathetic, 
effects. It increases (decreases) the electric field in a particular hemisphere, while 
simultaneously decreasing (increasing) the drift distance, thereby affecting the drift 
time. This effect was seen as the oscillatory structure in the left panel of Figure IB. 51 
Fortunately, it is possible to correct for this effect. 
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Figure B.5. The time position of clusters at the outer radius of the FTPC 
as a function of hardware sector. The affect of the inner cathode offset is 
clearly seen as an oscillatory structure in both FTPCs (0-30, FTPC W, 
30-60, FTPC E) (Left panel). After corrections are applied, the magnitude 
of the oscillations is reduced (Right panel). 
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